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PREFACE 


The study described herein was performed at the request of the US Army 
Engineer District, Detroit, and was authorized on 5 April 1983. 

The study was conducted from May 1983 to September 1984. From May 1983 
to July 1983, the study was performed at the US Army Engineer Waterways Exper- 
iment Station (WES), Hydraulics Laboratory (HL), under the direction of 
Mr. H. B. Simmons, Chief, Mr. C. E. Chatham, Acting Chief, Wave Dynamics Divi- 
sion (WDD), and Mr. D. G. Outlaw, Acting Chief, Wave Processes Branch. From 
July 1983 to September 1984, the study was performed at WES in the Coastal 
Engineering Research Center (CERC), under the direction of Dr. R. W. Whalin, 
Chief, Dr. Fred Camfield, Acting Chief, Engineering Development Division 
(EDD), and Dr. Dennis R. Smith, Chief, Prototype Measurement and Analysis 
Branch. CPT George M. Horsham conducted the study and prepared this report. 

Commanders and Directors of WES during the performance of the study and 
preparation and publication of this report were COL Tilford C. Creel, CE, and 
COL Robert C. Lee, CE. Technical Director was Mr. Fred R. Brown. 
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CONVERSION FACTORS, NON-SI TO SI (METRIC) 
UNITS OF MEASUREMENT 


Non-SI units of measurement used in this report can be converted to SI (metric) 


units as follows: 


Multiply By To Obtain 
feet 0.3048 metres 
miles (US statute) 1.609347 kilometres 
pounds (mass) 0.4535924 kilograms 
slugs (mass) per cubic foot 515.3788 kilograms per cubic metre 
3 
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WAVE CLIMATOLOGY STUDY FOR LUDINGTON HARBOR, MICHIGAN 


PART I: INTRODUCTION 


1. At the request of the US Army Engineer District, Detroit, the 
Coastal Engineering Research Center (CERC) of the US Army Engineer Waterways 
Experiment Station (WES) submitted a proposal for conducting a wave climatology 
study of Ludington Harbor, Michigan. The proposal was accepted, and funds for 
the study were authorized on 5 April 1983. 

2. The primary objective of the study was to measure wave conditions in 
Lake Michigan in the vicinity of the harbor entrance and, simultaneously, in 
the Ludington Harbor channel. The acquired wave data were then used to char- 
acterize the relationship between wave conditions in the inner harbor with 
those occurring in the lake. 

3. Ludington Harbor, Michigan, is located on the eastern shore of Lake 
Michigan, about 153 miles* northeast of Chicago, Illinois, and 60 miles north 
of Muskegon, Michigan (Figure 1). The harbor is configured as an outer basin 
formed by two shore-connected arrowhead breakwaters and an inner channel which 
connects the outer basin with the northern end of Pere Marquette Lake where 
the lake ferry berthing facilities are located (Figure 2). 

4. The scope of work of the study included design, deployment, 
maintenance, and recovery of two wave measurement systems, reduction and 
analysis of acquired wave data, correlation of observed wind and wave con- 


ditions, and evaluation of results. 


a 


* A table of factors for converting non-SI units of measurement to SI 
(metric) units is presented on page 3. 
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Figure 2. Wave gage locations 
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PART II: WAVE MEASUREMENT SYSTEM 


5. Each of the two wave measurement systems employed in the study 
consisted of a tripod instrumentation platform, wave gage, acoustic location 
device, steel mooring cable and float, and a lighted witness buoy. Figure 3 
is a schematic depiction of the system. The wave gages are self-contained, 
internal recording pressure sensing instruments. The gages were configured to 
acquire 1,024 pressure samples N at a rate of 1 Hz, resulting in record 
lengths of 1,024 sec TI , or about 17 min. A wave record was acquired every 
2 hr. 

6. The wave measurement systems were initially deployed on 26 May 1983 
with the assistance of Grand Haven Area Office, CE, personnel. Gage 9 was 
deployed in Lake Michigan; Gage 7 was deployed at the east end of Ludington 
Harbor Channel. Site selection for placement of the gages was made after 
reviewing information presented in the physical model study of Ludington 
Harbor.* The gages were serviced at approximately 60-day intervals. During 
the October 1983 servicing, Gage 9 was found to be defective and was replaced 
by Gage 67. The gages were retrieved from the measurement sites on 
10 December 1983. 

7. Table 1 provides a wave data statistical analysis summary, and 
Tables 2~4 provide a wave data record summary for Gages 7 and 9. Significant 
wave height versus time data for Gages 7, 9, and 67 are presented in Plates 


1-6. 


* L.G. Crosby and C. E. Chatham, Jr. 1975 (Sep). "Design of Entrance Chan- 
nel Improvements for Ludington Harbor, Michigan; Hydraulic Model Investiga~ 
tion," Technical Report H-75-14, US Army Engineer Waterways Experiment 
Station, Vicksburg, Mississippi. 
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PART III: DATA PROCESSING AND ANALYSIS 


8. The raw data were examined and edited for bad data points. Those 
wave records which contained greater than a selected small percentage of bad 
data points were discarded. The edited wave records then were spectrally ana- 
lyzed via a Fast Fourier Transform algorithm. The appropriate depth response 
factor was applied to each frequency band, and pertinent statistics were com- 
puted. The computed statistics include significant wave height H » peak 
period Le » and spectral density plots for selected wave records (Appendix B). 
The spectral plot scaling factor is constant to allow direct intercomparison 
of results. The plots are useful for intercomparison of wave conditions in 
the lake and channel at corresponding times. Appendix A provides details of 
the wave data processing procedure. 

9. Wind direction, wind speed, and significant wave height measured at 
both locations were correlated as functions of time during periods of signifi- 
cant wave activity. The correlation procedure consisted of plotting wind- 
speed, wind direction, and significant wave heights as ordinate values and 


times as the abscissa (Appendix B). 
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PART IV: CONCLUSIONS 


. 10. Evaluation of the data presented herein led to the following 
- conclusions. 

ll. Significant wave heights in Lake Michigan during the period of 
observation ranged from less than 1.0 ft to a maximum of 8.8 ft. Signifi- 
cant wave heights in Ludington Harbor Channel ranged from less than 1.0 ft to 
7 a maximum of 1.3 ft. The greatest peak period observed in Lake Michigan during 
the period of observation was 8.8 sec; the greatest peak period observed in 
the Ludington Harbor Channel ranged from about 6.2 sec to about 6.5 sec.* 

12. Inspection of peak periods occurring simultaneously in the lake and 
harbor channel indicate a very close relationship between the periods up to a 
maximum of about 6 sec. During times when the peak period in Lake Michigan 
exceeds 6 sec, the period in the harbor channel is still about 6 sec. It 
. appears the geometrical configuration of Ludington Harbor does not permit 
waves of periods greater than about 6 sec irrespective of the periods of the 
waves entering the harbor from Lake Michigan. 
fp 13. Comparison of the wind and wave records indicates that periods of 
z greatest wave activity in the harbor channel coincide with periods when the 

winds are from almost due west. It is during these periods that the peak of 
the wave spectra obtained in the harbor channel closely match the peak of the 
P spectra obtained in Lake Michigan; thus it is reasonable to assume that the 
waves observed in the channel are propagating in from Lake Michigan. However, 
” when the wind is from directions other than very close to due west, there is 
, almost no wave activity observed in the harbor channel even though the signif- 
y icant wave height in Lake Michigan just outside the harbor exceeds 4 ft. Con~ ‘ 
u versely, when the wind is from almost due west, the significant wave height in 
the harbor channel approaches the observed maximum of 1.3 ft even though the 


a significant wave height in Lake Michigan at the time is only about 2 ft. 


r. * The periods of 227.556 shown in Tables 2, 3, and 4 are a consequence of 
, spectrally analyzing wave records obtained during calm conditions and are 
ignored for the purpose of this discussion. 
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APPENDIX A: WAVE DATA PROCESSING PROCEDURE 


Subsurface Pressure 


1. Subsurface pressure under a wave is the summation of two contribut- 


ing components, dynamic and static pressure, and is quantitatively defined as ‘ 


7 cosh [2n(z + d)/L] H 2nx 21t 
P o8 cosh (2nd/L) 7 aerial Mah F T | ~ 082 * Pa (Al) 
where 
p' = total or absolute pressure 


pe = w/g = mass density of water (for saltwater, p = 2.0 1b- 


see te = 2.0 Stuga/eee ; for freshwater, p = 1.94 sluas/tt>) 


d = water depth 

H = wave height 

L = wave length 

T = wave period 
Pe atmospheric pressure 


The first term of Equation Al represents a dynamic component due to accelera- 
tion, while the second term is the static component of pressure. For conven- 


ience, the pressure is usually taken as the gage pressure defined as 


Equation A2 can be written as 
roan Sma BU “ 
2. It is often necessary to determine the height of surface waves based : 


on subsurface measurements of pressure. For this purpose it is convenient to 


rewrite Equation A3 as 


Al 


CNN NS INTE Tat tec anata teenie tard at Satan ati oa aan ae UA SET 


n® Mp one) (A4) 
Zz 


where 
N = correction factor equal to unity if the linear theory applies 
z = depth below the still-water level of the pressure gage 


K, = pressure response factor 


Several empirical studies have found N_ to be a function of period, depth, 
wave amplitude, and other factors. In general, N decreases with decreasing 
period, being greater than 1.0 for long-period waves and less than 1.0 for 
hort =pecioa waves. 


3. The pressure response factor K is therefore written as 


z cosh (kD) 


where 
k = local wave number 
= local water depth 
D = height of the sensor above the bottom 


4. The spectral energy of the sea surface is thus related to the pres~ 


sure spectrum by the following equation: 


cosh (kd) 


2 
E,(£) 7 Exc BE, (£) (A6) 


where subscripts s and p_ refer to the surface and pressure spectra, 


respectively. 


5. The significant wave height H, is quantitatively defined as 


l 1/2 


2At 


Hy = 4 f E,(£)df 
1/T 


The processed data were smoothed by band averaging eight frequency components 


of the raw periodogram. The final bandwidth resolution Be of the wave 


A2 
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spectra therefore is B, =>, = 0.0078 . 


Statistical Analysis 


6. The wave record data were analyzed to obtain statistical parameters 
as functions of wave frequency. The statistical analysis output parameters 
were 

th 
Band number - the n band of the wave frequency spectrum. 


a 
b. Wave frequency - the central frequency of each band in the 
frequency spectrum. 


Wave period - the period corresponding to the wave frequency. 


c 
d. Valid observations - the term "valid” means that the wave record 
contained 20 or fewer bad pressure samples. 


e. Percent of valid records - the quantity of valid observations 
compared to the total quantity of wave records. 


f. Probable maximum wave height - the probable height of the 
highest wave predicted by the Rayleigh Wave Height probability 
distribution, a height ratio of 1.87. 


A3 


APPENDIX B: WAVE DATA SUMMARY 


1. This Appendix presents a summary of wave data. Pages B2, B25, B51, 
and B62 present wind speed, wind direction, and significant wave height data 
plots for (a) 13 October 1983 - 15 October 1983, (b) 20 November 1983 - 

22 November 1983, (c) 24 November 1983 - 26 November 1983, and (d) 29 November 


1983 - 1 December 1983, respectively. Corresponding spectral plots are pre- 


sented for time periods (a), (b), (c), and (d) on the following pages: 
(a) pages B3-B24, (b) pages B26-B50, (c) pages B52-B6l, and (d) pages B63-B82. 
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